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Watanabe and Hu Reply: In a recent Letter 
Watanabe, Yukawa, Ito, and Hu (WYIH) proposed the 
idea of superscaling for percolation on rectangular do- 
mains and used numerical data near the effective critical 
density to support such idea. In |^, Pruessner and 
Moloney (PM) pointed out some problems in Most 
of these problems are caused by the differences in using 

or the bulk critical density pc in the analysis. Figure 
2 of Q shows clearly that at p = pc the existence prob- 
abilities Ep for different aspect ratios R are different Q- 
Thus it is natural that we consider a horizontal line at 
Ep = 0.5 in Fig. 2(a) to define the effective critical den- 
sity p'c for curves with different aspect ratios R and usiM 
data near to study superscaling shown in Fig. 4 of P . 
Piease note that at p(,, the correlation length is of the 
order of vertical linear dimension L. By considering Ep 
near rather than near pc, we will not have the prob- 
lems mentioned at the second and the third paragraphs 
of Ø- To avoid the confusion, in the abstract and the 
first paragraph of Q we should replace “deviation from 
the critical point e” by “deviation from the effective crit- 
ical point p'c, e,”; we should also replace the term “the 
critical point” below Eq. (1) and “pc” of Eq. (2) by “the 
effective critical point p(,” and , respectively. 

The value of b in Eq. (10) of Q was determined by PM 
to be 3/8(= yt/2) based on P{L,R,e) = P{RL,R~^,e), 
while WYIH obtained b = 0.05. This discrepancy is 
caused by differences in definitions of percolating clus- 
ters. PM used only the largest cluster to define the per- 
colation probability P and thus obtained P{L,R,e) = 
P{RL, R~^ ,e). WYIH used all clusters extending from 
the top to the bottom of the system to define P and have 
P{L,R,e) > P{RL,R-\e) for R > 1. Resides WYIH 
used the effective critical density p(, in Fig. 7 of [y to de- 
termine b = 0.05, while PM did not used p(,. To avoid the 
confusion, we should also use the effective critical density 
rather than the critical density to define t in Eq. (10) of 
P as in the case of Fig. 7. 

Since a horizontal line at Ep = Ep{L, 0,1) in Fig. 2 
of 0 has been used to defined p(., Fig. 2 indicates that 
p(, depends on L and R and limi^oo Ep[L, p(, — pc, R) — 
Ep{L,0,R) ^ 0 for R ^ 1. Thus PM correctly pointed 
out that the difference between using p(. and pc does not 
vanish in the limit L — > oo for R ^ 1. For this reason, we 
should replace “correct for finite-size effects” above Eq. 
(9) in Q by “study superscaling behavior” . 

The meaning of the nonlinear scaling variable ques- 
tioned by PM can be understood from the mapping func- 
tion /(p) defined by Ep{L,p,R) = Ep{L, f{p),l), with 
/ determined numerically and plotted in Fig. 1 which 
shows nonlinear dependence on p for i? > 1 and is con- 


sistent with Fig. 6 and Table I of Q. Finite-size scal- 
ing implies Ep{L,f{p),l) ~ F{{f{p) - pc)Ly^); com- 
bining this with the expansion of /(p) at p(. up to the 
linear term with {df/dp)p'^ ~ i?“ gives Ep{L,p,R) ^ 
E{{p — p'^Ly*R°-), which shows superscaling of Ep with 
scaling variable x = (p — Pc)Ly^R°‘ and p(, (instead of pc) 
enters x naturally. We find that ln(9//i9p)p^ increases 



P 

FIG. 1: The mapping functions / in Eq. (1), which are non- 
linear for R > 1. The inset shows how to obtain / from Ep. 
The solid circles are the effective critical points p(. 


linearly with In i? with a well defined exponent a = 0.14 
0 which is consistent with that obtained in P|. Similar 
calculations of the mapping functions of the percolation 
probability give the exponent b — 0.05 0 which is also 
consistent with that in P|. 

Scaling can be considered as “data collapse” for critical 
systems with different sizes 0 . We found very good data 
collapse for systems with different aspect ratios in P| . It 
is appropriate to call such results “superscaling” . 
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